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We consider a model for bilayer superconductors where 
an interlayer pairing amplitude Ax co-exists with intralayer 
pairing Ay of d x 2_ y 2 symmetry. This model is motivated by 
a recent photoemission experiment reporting the splitting of 
the nodes of the energy gap. In addition to offering a natural 
explanation of this observation, the model has a number of 
new experimental consequences. We find that the new state 
is accompanied by a spontaneous breaking of the tetragonal 
symmetry. We also find that the out-of-phase oscillation of 
Ax and Ay gives rise to a new Raman active mode. The phase 
of Ax may also become imaginary, leading to a state which 
breaks time reversal symmetry, which may have important 
implications for tunnelling experiments. 

PACS numbers: 74.65. +n, 74.50. +r 

In the past two years, there has been mounting evi- 
dence for nodes in the gap in oxide superconductors . 
Recently, a new class of experiments which are phase sen- 
sitive have provided strong evidence for a sign change of 
the order parameter as a function of angle consis- 
tent with d symmetry. At the same time, there are sev- 
eral experiments which appear to be inconsistent with 
the simple d x 2_ y 2 state. For example, Sun et al Q, have 
made tunnel junctions between Y Ba^Cu^O^-s (YBCO) 
and Pb and obtained Josephson current along the c— axis. 
Secondly, a recent report of fractional vortices tied to 
grain boundary interfaces |6) has been interpreted as re- 
quiring a time reversal symmetry breaking state near the 
interface 0. Finally, Ding et al J§| measured the en- 
ergy gap as a function of angle usng angular resolved 
photoemission in Bi^Sr^C 'aCuiO% (Bi-2212) and claim 
that the energy gap does not vanish along (n, tt) as ex- 
pected for the d x 2_ y 2 state. Instead, the node is split into 
two, lying approximately 10° from the 45° line. Another 
group [^| also reported a finite energy gap in the (tt, tt) 
direction and studied its temperature dependence. Taken 
together, these experiments suggest that the energy gap 
structure may be more complicated than d x 2_ y 2. In this 
paper we show that a modification of the d x 2_ y 2 state 
due to interlayer pairing provides a scenario whereby the 
experiments mentioned above can be explained. Our pic- 
ture also leads to a number of predictions which can be 
tested experimentally. 

Band structure calculations have yielded a surprisingly 
large interlayer hopping term between the bilayers, given 
by t±(k)c\ rr (k)c2 <7 (k) where 1 and 2 refer to the layers. 
Chakravarty et al |10f| have proposed the parametriza- 



tion tx(fc) = ^^(cosA^a — cosfc y a) 2 with t° ± w 0.24 eV. 
Andersen et al [0 recently showed that the unusual k de- 
pendence originates from hopping via the Cu 4s orbitals, 
and that it leads to an interlayer exchange energy of order 
Jx ~ 20 meV, consistent with the lower bound of 8 meV 
given by neutron scattering fl2|| . It has been proposed 
jl3],[l4[ that the relatively large exchange term, when en- 
hanced by intralayer antiferromagnetic correlations, is 
responsible for the spin gap phase found in bilayer sys- 
tems [^5). Furthermore, Ubbens and Lee [[l3| have argued 
that in the superconducting phase, a pairing amplitude 
A12 = (ci^(r)c2i(r)) appears in addition to the intra- 
plane order parameter An (if) = (a^(r)ci[(r+ if)), i = 1,2 
which is of d symmetry. This leads to a quasi-particle 
dispersion relation E±(k) = (£| + |A|| ± Ax| 2 ) 5 where 
6c = (-k — ^, Ay (k) and Ax. are proportional to An = A22 
and A12 respectively. Since the nodes of this state are 
given by the zeros of |A|| ± Aj_|, if A|| has d x 2_ y 2 symme- 
try, the nodes are shifted from the 45° direction and split 
into two nodes Jl3| . This provides a natural explanation 
of the observation by Ding et al 

Encouraged by the experiment, we decided to re- 
examine the interlayer pairing model. Ubbens and Lee 
]l3| included Jj_ but not the tj_ term in their considera- 
tion. In this paper we add both the t± and Jj_ terms to 
the standard t — J model. Wheras Ubbens and Lee at- 
tempted to justify the appearance of A12 microscopically, 
in this paper we take a more phenomenological approach 
and assume the co-existence of Ax and An. This is mo- 
tivated by experiment: as far as we know, the present 
scenario is the only one which is consistent with both the 
node splitting [|| and the sign change of the order param- 
eter as k varies from to 7r/2, as required by the corner 
SQUID experiment The purpose of this work is to 
explore the consequences of the assumed interlayer pair- 
ing amplitude, so that further experiments can confirm 
or falsify this picture. 

We begin by treating the bilayer t — J model (includ- 
ing t± and J±) using the slave boson mean field method, 
which is equivalent to the Gutzwiller approximation. The 
tight binding band e(k) = — 2t(cos k x a + cos k y a)) is split 
into bonding and anti-bonding bands g± — (ci ± c-i )/ \/2 
with dispersion e±(k) = e(k) ± t±(k), where t is of or- 
der J and t±(k) ~ Xot±(k) where Xq is of order x, the 
doping concentration. Thus, in the mean field theory the 
effective interlayer hopping is reduced by x, simply be- 
cause in the strongly correlated metals, the electron, on 
average, must find a vacancy to hop. 
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The resulting Fermi surface is shown in Fig 1 for x — 
0.15. As emphasized by Anderson |l6|, in the normal 
state, coherent hopping is really not possible between 
the layers. We expect the bonding-antibonding splitting 
to be smeared out, but a region of low lying excitations 
may exist in the k space between the two Fermi surfaces. 
Fig. 1 bears a striking resemblance to the photoemission 
data of Dessau et al |l7]] . 

In the superconducting state, coherent hopping be- 
tween the planes indeed occurs, and we should take 
the band splitting seriously. With the basis sets 



( C l,feT' C i ^ki ,C 2,kV C 2-kl J 

takes the form 



the mean field Hamiltonian 



e(k)-n A,|(fc) t±(k) A ± 

A|(fc) -e(k)+fi A* ± -t x (k) 

t±(k) Aj_ e(fe)-/" A||(fc) 

-h(k) A|(*) -e(k)+fi 



A* 



(1) 



It is easily seen that the bonding anti-bonding bands g± 
block diagonalize this matrix, so that the g± bands are 
separately paired by A±(fc) = A||(fc) ± Aj_, resulting in 
the quasi-particlc spectrum 



E±(k) = (e±(k) + \A ll (k)±A_ 



(2) 



where £±(fc) = e±(fc) — We can choose A|| to be real 
and positive. If Aj_ is also real and positive, the nodes 
are split as before, but the split nodes are associated with 
the ± bands separately, as indicated in Fig. 1. It is clear 
from Fig. 1 that for real Aj_, the onset of interlayer 
pairing implies a spontaneous breaking of the tetragonal 
symmetry of the model. This is a consequence of the 
inclusion of t± . The two degenerate states of the broken 
symmetry correspond to Aj_ being positive or negative. 
We shall refer to these degenerate states as d± s. Due 
to the vanishing of t±(k) along (-7T, 7r), this assymmetry 
is difficult to resolve near the nodes. However, in a given 
domain, the g+ band has a different gap along MY than 
MY and the resulting asymmetry in the electronic state 
is expected to couple linearly to the orthorhomic strain 
e=(a-b)/(a + b). 

There is one additional complication to this discus- 
sion, in that we should consider the possibility that Aj_ 
is purely imaginary. We shall refer to this state as d + is, 
which has a minimum gap of |A_l|. Indeed, a mean 
field calculation shows that this state is lower in en- 
ergy than d + s. The Ginzburg-Landau free energy is 
given by F = Fq + F e , where for simplicity we have fixed 
An = A 22 = A,,, 



= ay | A,, 



2 + h I A,, | 4 
di|A|,| 2 |A^ 



IA; 



&±|A, 



d 2 (A|Ai* +c.c.) (3) 



where F e describes the linear coupling to the strain dis- 
cussed earlier. jl9| In mean field theory, we found that 
di = 4d 2 > 0. If we write Aj_ = |AJe^, the d 2 term 
is proportional to cos 2cf> with a positive coefficient, and 
is minimized by (f> — ^. The F e term, on the other 
hand, is proportional to cos</>, so that </> = or ir (the 
d ± s state) is stabilized in systems with pre-existing or- 
thorhomic distortions such as YBCO. In tetragonal sys- 
tems, we add an elastic energy term Fe = \k£ 2 to F 
and minimize with respect to e. This produces a term 
v 2/„\ia |2| a .. 1 2 „„c2 a wmc h opposes the d 2 term. 



-2(a 2 /K)|A ± | 2 |A,|| 2 cos 2 1 



F e = ae(A||A^ + c.c.) 



(4) 



The d + s state is stabilized provided ^l. > 2d 2 ■ The ob- 
servation of split nodes in the nominally tetragonal Bi- 
2122 presumably means that this condition is satisfied 
and d + s is stabilized. Thus we predict that, for bilayer 
systems with tetragonal symmetry, a superconducting 
state with split nodes is accompanied by a spontaneous 
breaking of the tetragonal symmetry. The Bi — 2122 
compound is not truly tetragonal due to the superlattice 
modulation in the (7r, 7r) direction, but the a and b lat- 
tice constants are predicted to become unequal at low 
temperature. The best test of the prediction is proba- 
bly in tetragonal materials such as the bilayer mercury 
compounds. Even in YBCO, we expect an additional 
distortion below T c . There is evidence for this in the 
literature |l8|. The coupling to lattice distortions also 
lead us to expect anomalies in the transverse ultra sound 
velocity below T c Q. 

As mentioned earlier, the d + s state explains the pho- 
toemission data of ref. 8. Indeed, the lower branch min 
(A + , A_) shown in the insert of Fig. 1 can be fitted to 
the experimental data with A||(fc = 0) ~ 30 meV and 
A± K 5 to 9 meV pb]| . In ref. 8, only data associated 
with the g_ band (dashed line in Fig. 1) was shown. The 
question arises as to what happened to the upper branch 
of the insert in Fig. 1. Two possibilities need to be exam- 
ined: the photoemission may be from a single domain of 
d + s state or from multiple domains of d+s and d— s. In 
the first case we expect a single peak with different gaps 
(A||(fc = 0) ± A_l) at the M - Y and M -Y crossings. 
In the second case the peak should be a superposition of 
two peaks with a splitting of 2A_l. The data is consistent 
with a splitting of 10 meV or less p(| which may require 
a slight angular dependence of A^(fc) so that it is smaller 
along r — M than along T — Y. At present, data is not 
available which covers both T — M and T — M quadrants 
in the same sample, so that these two possibilities cannot 
be distinguished. 

The existence of several order parameters should lead 
to new collective modes, which are amplitude and phase 
modes associated with Ajj = \Aij\e l ^ ij ,i, j = 1,2. The 
amplitude modes are expected to be high in frequency 
and damped and we shall focus on the phase modes only. 
There are three modes corresponding to the three phase 
degrees of freedom. 
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(1) The in-phase oscillation of An,A22,Ai2 is the 
Bogoliubov-Anderson mode which is coupled to total 
charge density and pushed up to the plasma frequency. 

(2) The out-of-phase mode <fiu — <j>22- Recently this 
mode was discussed |2lJ in the context of the interlayer 
tunnelling model of Chakravarty et al [JlO|. However, 
as pointed out by Wu and Griffin, the existence of this 
mode is a general feature of two layers coupled by Joseph- 
son tunnelling and was described by Leggett p2] ], who 
noted that the usual phase-number commutation rela- 
tion [4>u,Ni] — 2i implies a coupling of </>n — ^22 to 
Ni — N2, i.e., to charge transfer between the planes. The 
mode frequency luj can be computed by combining the 
Josephson relations <f> — 2 eV and h\2 = Ej sin <j> with 
eV — 7ni2, where 7 = C^ 1 + dfx/dn. n±2 is the inter- 
layer density oscillation, C is the capacitance per area, 
and Ej « |fj_| 2 A^ 2 (0)|A|| \a 2 is the Josephson energy per 
area between the planes. We find wj = (2Ejj)i which 
is just the Josephson plasma frequency. If we ignore 
the capacitance term and approximate 7 = dfi/dn by 
iV(O) -1 m Ja 2 , wc find that ujj « xt±(A\\/J)i which 
is quite high, so that this mode is probably strongly 
damped. 

(3) Finally, the mode which is of greatest interest to us 
is the out-of-phase oscillation between c/>x2 and <j)n = 4>2i- 
Since the free energy given by Eq. (3) depends weakly on 
this phase difference, we expect this mode to be low-lying 
in frequency. Furthermore, this phase mode leads to a re- 
arrangement of the quasi-particle spectrum in the plane, 
since a change in </> 12 from to ^ interpolates between 
the d + s and d + is states. This should couple to charge 
oscillation within each plane, but with the total charge 
conserved within each layer. Thus, we expect this mode 
to be Raman active. To confirm this qualitative picture, 
we have carried out a collective mode calculation in a 
model where we introduce phcnomcnological attractive 
coupling constants to stabilize A|| and Aj_- We treat 
this model in mean field theory and carry out an expan- 
sion around the mean field in a standard way For 
simplicity we have set e = 0. We find that the globally 
stable states correspond to cither A|| or Aj_ = 0, while 
the d+is state lies slightly higher in energy but is locally 
stable. The global energy balance depends on details of 
the interaction and the presence of e. For example, a 
sufficiently large e will stabilize the d + s state because 
of the linear coupling in Eq. (4). Since our goal in the 
part of the analysis is to understand the collective mode 
rather than the microscopic origin of the energy gaps, we 
proceed to expand about the locally stable d + is state. 
We find that the Josephson mode </>n — $22 decouples as 
expected. Setting <pu = ^22, we find that <fiii,<fii2 and 
P2 are coupled, where P2 is the I = 2 component of the 
density fluctuation in the plane. The 3x3 matrix can 
be diagonalized. Details of this calculation will be given 
elsewhere, but the main result is an estimate of the out- 



of-phase mode frequency, which turns out to be of order 
|A_l|- For |Aj_| <C I A || I, its damping should be small. 
Furthermore, its coupling to P2 implies that it is Raman 
active in the B\ g symmetry [^J . The detection of this 
mode will be a strong confirmation of the existence of 
A x . 

Next we examine what effects the more complicated or- 
der parameter structure have for tunnelling experiments. 
In ref. 4, tunneling is between oxide superconductors 
across a grain boundary. A reasonable tunnelling Hamil- 
tonian may be of the form H t = to c iL°iR + c - c -> 

i=l,2 

which will lead to a Josephson coupling of the form 

Hj = Tj(A* +L A +R + A*_ l A-r + c.c.) (5) 

The two terms in Eq. 5 will generate two Josephson 
currents with different dependences on the grain bound- 
ary orientation f22j ] . However, at present the half integer 
flux experiment in specially designed rings Q is done on 
twinned samples, so that the d ± s domains are difficult 
to distinguish from the pure d x i_ y i state. 

In the experiment which measures the critical current 
of a SQUID loop ||, the tunnel junctions are between 
the oxide superconductor and a conventional supercon- 
ductor, and the tunnelling current is in the a — b plane. 
In this case, it may be more reasonable to assume a tun- 
nelling Hamiltonian of the form H' T = tog+cji+c.c. where 
cr destroys an electron in the conventional superconduc- 
tor. In this case, only the symmetric order parameter 
A + forms the Josephson coupling. Thus, in principle, if 
a single domain sample can be made with variable junc- 
tion orientations, the shift of the node shown by the solid 
line in the insert of Fig. 1 can be detected. 

We briefly discuss the experiment of ref. 5, where a 
Josephson current is observed along the c axis between 
YBCO and a conventional superconductor. The d + s 
state by itself cannot explain why the critical current is 
relatively insensitive to whether the YBCO is untwinncd 
or highly twinned. This is because the coupling to the 
lattice strain e (shown in Eq. (4)) will lock the d + s state 
to one set of twins and the d— s states to the other. Since 
the d order parameter is insensitive to twinning (other- 
wise the phase sensitive experiment of references 3 and 
4 will not work in twinned samples), the Josephson cur- 
rent will have opposite signs on the two sets of twins and 
tend to cancel. One possible way out of this dilemma 
is to postulate the existence of regions where the d + is 
becomes stabilized. This is plausible due to the small en- 
ergy differences expected and can happen near the twin 
boundaries. With this assumption a net Josephson cur- 
rent may be obtained. We defer a detailed discussion of 
this possibility to a later publication |24|] , 

The trapped fluxes on grain boundaries between 
YBCO films with different orientations are measured and 
interpreted as being due to the appearance of fractional 
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charged vortices ||. Such vortices require the existence 
of a time-reversal symmetry breaking state near the in- 
terface 0. Recently, Kuboki and Sigrist ]2j| offered an 
explanation of this state as being due to an admixture 
with a proximity induced s component of the order pa- 
rameter. In our picture, the possible existence of the 
complex d+is state near an interface suggests an alterna- 
tive origin of the time-reversal symmetry breaking state. 
The two alternatives can be distinguished by searching 
for similar effects in single layer materials. 

Finally, we address the issue of the possibility of a sec- 
ond phase transition below T c . Usually, a Landau theory 
of the form given by Eq. (3) with two order parameters 
implies two phase transitions, because the coefficients a± 
and an change signs at different temperatures. In the 
underdoped case, the transition is complicated by strong 
fluctuation effects, so that, in some sense, the spin gap 
transition may be considered the first transition, and the 
superconducting T c the second transition. This raises the 
interesting possibility that the situation may be reversed 
in optimally doped or overdoped cases, which would im- 
ply that the transition at T c is to a pure d x i_ y i order 
paramter, followed at some lower temperature by the on- 
set of interlayer pairing (Aj_ 7^ 0). Recently, a second 
transition at is reported in T2 measurements in fully 
oxygenated samples of YBCO po| . Obviously, it will be 
interesting to look for lattice distortions near this tem- 
perature. Another interesting possibility is to study the 
fate of the fractional vortices at temperatures above 30K . 
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Figure 1. Mean field calculation of the Fermi sur- 
face including hopping between bilayers for hole doping 
x = 0.15. The parameters t/J = 3,tf } /J = 2.4 and 
J±/J = 0.2 are used. The solid and dashed lines are 
the Fermi surfaces for the bonding (g+) and antibonding 
(<7_) bands. Solid and dashed arrows indicate the ap- 
proximate location of the nodes associated with A + and 
A_ in the d + s state. The insert shows a schematic pic- 
ture of the angular dependence (0 = is along the T — M 
axis) of the gap functions |A+| (solid) and |A_| (dashed) 
for the d + s state (Aj_ real and positive). Note that A± 
are associated with the solid and dashed Fermi surfaces 
respectively. For the d — s state, |A+| is given by the 
dashed line. 
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